ANISOTROPIC STEP STIFFNESS FROM A KINETIC MODEL OF 

EPITAXIAL GROWTH 



DIONISIOS MARGETIS* AND RUSSEL E. CAFLISCHt 

Abstract. Starting from a detailed model for tiie kinetics of a step edge or island boundary, we 
derive a Gibbs-Thomson type formula and the associated step stiffness as a function of the step edge 
orientation angle, 9. Basic ingredients of the model are: (i) the diffusion of point defects ("adatoms") 
on terraces and along step edges; (ii) the convection of kinks along step edges; and (iii) constitutive 
laws that relate adatom fluxes, sources for kinks, and the kink velocity with densities via a mean- 
field approach. This model has a kinetic (nonequilibrium) steady-state solution that corresponds 
to epitaxial growth through step flow. The step stiffness, /3{6), is determined via perturbations of 
the kinetic steady state for small edge Peclet number, P, which is the ratio of the deposition to the 
diffusive flux along a step edge. In particular, /3 is found to satisfy = O{0~^) for 0{P^^'^) < 9 <C f , 
which is in agreement with independent, equilibrium-based calculations. 
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1. Introduction. The design and fabrication of novel small devices require the 
synergy of experiment, mathematical modeling and numerical simulation. In epitax- 
ial growth, crystal surface features such as thin films, which are building blocks of 
solid-state devices, are grown on a substrate by material deposition from above. De- 
spite continued progress, the modeling and simulation of epitaxial phenomena remains 
challenging because it involves reconciling a wide range of length and time scales. 

An elementary process on solid surfaces is the hopping of atoms in the presence of 
line defects ("steps") of atomic height [T31 HSJ |31] : atoms hop on terraces, and attach 
to and detach from step edges (or island boundaries). Burton, Cabrera and Frank 
(BCF) [6 first described each step edge as a boundary moving by mass conservation 
of point defects ("adatoms") which diffuse on terraces. In the BCF theory, the step 
motion occurs near thermodynamic equilibrium. Subsequent theories have accounted 
for far-from-equilibrium processes; for a review see section [2j 

The macroscale behavior of crystal surfaces is described by use of effective mate- 
rial parameters such as the step stiffness, /3 [27]. In principle, f3 depends on the step 
edge orientation angle, 0, and is viewed as a quantitative measure of step edge fluctu- 
ations [TJ [3H] • Generally, effective step parameters such as /3 originate from atomistic 
processes to which inputs are hopping rates for atoms; in practice, however, the pa- 
rameters are often provided by phenomenology. For example, the dependence of f3 on 
9 is usually speculated by invoking the underlying crystal symmetry [H [191 [Ml |35] . 

In this article we analyze a kinetic model for out-of-equilibrium processes [3 H] in 
order to: (i) derive a Gibbs-Thomson (GT) type formula, which relates the adatom 
flux normal to a curved step edge and the step edge curvature [IT] |22]; and (ii) 
determine the step stiffness (3, which enters the GT relation, as a function of 9. For 
this purpose, we apply perturbations of the kinetic (nonequilibrium) steady state of 
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the model for small Pcclet number P, which is the ratio of the material deposition 
flux to the diffusive flux along a step edge, i.e. 

(1.1) P = i2a^f)/DE , 

in which a is an atomic length, / is a characteristic size for the flux / normal to the 
boundary from each side, and De is the coefficient for diffusion along the boundary. 



A factor of 2 is included in ( 1.1 ) since the flux is two-sided and the total flux is of size 
2/. For sufficiently small 6 and P, we flnd that the stiffness has a behavior similar to 
that predicted by equilibrium-based calculations [57] . 

For the boundary of a two-dimensional material region, a definition of /3 can arise 
from linear kinetics. In the setting of atom attachment-detachment at an edge, this 
theory states that the material flux, /, normal to the curved boundary is linear in 
the difference of the material density, p, at the boundary from a reference or "equi- 
librium" density, po. The GT formula connects po to the boundary curvature, k. For 
unit layer thickness and negligible step interactions [22], the normal ffux reads 

(1.2) f = DA{p-po), 

where Da is the diffusion coefficient for attachment and detachment, and pQ is deflned 

by 

(1.3) po = p,eT^ ^ pJl + j^k], \^K\<t:kBT. 



The last equation is referred to as the GT formula, in accord with standard thermo- 



dynamics [5J [TTlUll [SI] • In (1-3 1, p* is the equilibrium density near a straight step 
edge and fc^T is Boltzmann's energy (T is temperature); the condition \(3k\ <C ksT 
is satisfied in most experimental situations (43j . Equation ( [1.2| does not account for 
step permeability, by which terrace adatoms hop directly to adjacent terraces [^HTj . 
This process is discussed in section [2] 



For systems that are nearly in equilibrium, the exponent in (1.3 1 is derived by 
a thermodynamic driving force starting from the step line tension /3, the free energy 
per unit length of the boundary [18^ . The step stiffness f3 is related to /3 by [H [T5j [16] 

(1.4) P^P + Pee {Pe-^deP). 



Evidently, the knowledge of /? alone does not yield /3 uniquely: by (1.4 1, 
(1.5) P{e) = Ci cos + C2 sin + / M /3(i9) sin(6l - 1?) 







where Ci and C2 are in principle arbitrary constants. 

The parameters /3 and /3 are important in the modeling and numerical simulation 
of epitaxial phenomena. In thermodynamic equilibrium, the angular dependence of 
the step line tension, P{6), determines the equilibrium (two-dimensional) shape of step 
edges or islands, e.g. the macroscopic ffat parts ("facets") of the step are found by 
minimizing the step hue energy through the Wulff construction [20 ] 129 ] [30 ] 140 ] l42l |45] . 
Near thermodynamic equilibrium, the step stiffness, P{9), controls the temporal de- 
cay of ffuctuations from equilibrium [TJ [25] . The significance of /3 was pointed out 
by de Gennes in the context of polymer physics almost forty years ago [101 [H]: the 
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energy of a polymer (or step edge) can be described by a kinetic energy term pro- 
portional to (3 ■ {dx / dy)'^ , i.e., the stiffness times a "velocity" squared where x and 
y are suitable space coordinates and y loosely corresponds to "time." Starting with 
a two-dimensional Ising model, Stasevich et al. [36, 37, 38] [31] carried out a direct 
derivation of P{6) and P{9) from an equihbrium perspective based on atomistic key 
energies. For most systems, however, there has been no standard theoretical method 
for determining (3(6) and (3(0). 

More generally, energetic principles such as a thermodynamic driving force are 
powerful as a means of describing the macroscopic effect of atomistic kinetics. The 
range of validity of energetic principles is not fully known and is an important unre- 
solved issue. We believe that energetic arguments should be valid for systems that 
are nearly in local equilibrium, where the relevant processes approximately satisfy 
detailed balance. For systems that are far from equilibrium, however, energetic prin- 
ciples may serve as a valuable qualitative guide, even if they are not quantitatively 
accurate. 

The kinetic and atomistic origin of a material parameter that plays the role of the 
step stiffness are the subject of this article. For a step edge or an island boundary on 
an epitaxial crystal surface, we use the detailed kinetic model formulated by Caflisch 
et al. [71 [H] and further developed by Balykov and Voigt [H [3] for the dynamics of 
the boundary. The basic ingredients are: (i) diffusion equations for adatom and edge- 
atom densities on terraces and along step edges; (ii) a convection equation for the 
kink density along step edges; and (iii) constitutive, algebraic laws for adatom fluxes, 
sources for kinks and the kink velocity by mean-field theory. This model admits a 
kinetic (nonequilibrium) steady state that allows for epitaxial growth via step flow. 
The model has been partly validated by kinetic Monte Carlo simulations [7]. 

The detailed step model described in [7l[8] and section 2.3 focuses on the kinetics 
of adatoms, edge-atoms and kinks at a step edge. As discussed by Kallunki and 
Krug , an edge-atom is energetically equivalent to two kinks. For example, the 
equilibrium density of kinks is proportional to exjp[—s/{kBT)] while the equilibrium 
density of edge-atoms is proportional to exp[—2e/{kBT)], in which e is defined as the 
kink energy in [23], or identified with —{kBT/2)log{DK / De) in [7 ; Dk and De are 
diffusion coefficients for kinks and edge-atoms. On the other hand, the kinetics in [7l|8] 
are different for edge-atoms and kinks, since edge-atoms can hop at rate De, while 
kinks move through detachment of atoms at rate Dk- This situation is consistent 
with the kinetics described in |23j , in which the De and Dk are proportional to 
exp[-Est/{kBT)] and exp[-Edet/{kBT)], respectively. 

We are aware that the mean-field laws applied here, although plausible and an- 
alytically tractable, pose a limitation: actual systems are characterized by atomic 
correlations, which can cause deviations from this mean-field approximation. In par- 
ticular, the validity of the mean-field assumption may be limited to orientation angles 
9 in some neighborhood of 6* = 0. Note also that the most interesting results of this 
analysis are for 9 near zero. Determination of the range of validity for this model 
is an important endeavor but beyond the scope of this paper. An extension of this 
model, which could improve its range of validity, would be to explicitly track the 
kinks in a step edge. This additional discreteness in the model would make the anal- 
ysis of step stiffness more difficult. Our analysis is a systematic study of predictions 
from the mean-field approach only, and the conclusions presented here are all derived 
within the context of this approach. On the other hand, our analysis is more detailed 
than previous treatments of step stiffness, since it is based on kinetics rather than a 



4 



D. MARGETIS AND R. E. CAFLISCH 



thermodynamic driving force. Moreover, the model includes atomistic information, 
through a density of adatoms, edge-atoms and kinks. 

For evolution near the kinetic steady state, we derive for the mass flux, /, a term 
analogous to the Gibbs- Thomson formula ( |1.3[) , and subsequently find the correspond- 
ing angular dependence of the step stiffness, /3{9). Our main assumptions are: (i) the 
motion of step edges or island boundaries is slower than the diffusion of adatoms 
and edge-atoms and the convection of kinks, which amounts to the "quasi-steady ap- 
proximation" ; (ii) the mean step edge radius of curvature, k^^, is large compared to 
other length scales including the step height, a; and (iii) the edge Peclet number, P, 
given by is sufficiently small, which signifies the usual regime for molecular beam 
epitaxy (MBE). To the best of our knowledge, the analysis in this paper offers the 
first kinetic derivation of a Gibbs-Thomson type relation and the step stiffness for all 
admissible values of the step edge orientation angle, 9. (This approach is distinctly 
different from the one in e.g. [33] where classical elasticity is invoked.) Our results for 
the stiffness are summarized in section [s] see (3.7l~ pJ9| . 

A principal result of our analysis is that (3 — 0{9~^) for 0{P^/^) < <C 1, which 
by (1.5 1 yields (3 = 0{9\ii9) for the step line tension. This result is in agreement with 
the independent analysis in |361I37I [351 139] . which makes use of equilibrium concepts. 
A detailed comparison of the two approaches is not addressed in our analysis. Our 
findings are expected to have significance for epitaxial islands, for example in predict- 
ing their facets, their roughness (e.g., fractal or smooth island boundaries) and their 
stability, as well as for the numerical simulation of epitaxial growth. More generally, 
our analysis can serve as a guide for kinetic derivations of the GT relation in other 
material systems. For example, it should be possible to derive the step stiffness for a 
step in local thermodynamic equilibrium within the context of the same model. This 
topic is discussed briefly in section [6] 

The present work extends an earlier analysis by Caflisch and Li |8j, which ad- 
dressed the stability of step edge models and the derivation of the GT relation. The 
analysis in [8], however, only determined the value of (3 along the high-symmetry ori- 
entation, 9 = Q. This restriction was due to a scaling regime used in [5] on the basis of 
mathematical rather than physical principles. In the present article we transcend the 
analytical limitations of [5] by applying perturbation theory guided by the physics of 
the step-edge evolution near the kinetic steady state. 

Our analysis also leads to formulas for kinetic rates in boundary conditions in- 
volving adatom fluxes. In particular, the attachment-detachment rates are derived as 
functions of the step edge orientation, and are shown to be different for up- and down- 
step edges. This asymmetry amounts to an Ehrlich-Schwoebel (ES) effect |TTl 13^ . 
due to geometric effects rather than a difference in energy barriers. In addition, if the 
terrace adatom densities are treated as input parameters, the adatom fluxes involve 
effective permeability rates, by which a fraction of adatoms directly hop to adjacent 
terraces (without attaching to or detaching from step edges) [HI UHl SI] • Our main 
results for the kinetic rates are described by (3.1l-(3.5l. 

In this article we do not address the effects of elasticity, which are due for instance 
to bulk stress. One reason is that elasticity requires a non-trivial modiflcation of the 
kinetic model that we use here. This task lies beyond our present scope. Another 
reason is that, in many physically interesting situations, the influence of elasticity 
may be described well via long-range step-step interactions that do not affect the step 
stiffness. The study of elastic effects is the subject of work in progress. 

The remainder of this article is organized as follows. In section [2] we review the 
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relevant island dynamics model and the concept of step stiffness: In section [2. 1| we in- 
troduce the step geometry; in section [2T2] we outhne elements of the BCF model, which 
highlight the GT formula; in section [2.3| we describe the previous kinetic, nonequi- 
librium step-edge model [TJISj, which is slightly revised here; and in section 2.4 we 
outline our program for the stiffness, based on the perturbed kinetic steady state for 
small step edge curvature, k. In section [3] we provide a summary of our main results. 
In secti on [4| we derive analytic formulas pertaining to the kinetic steady state: In 
section |4.2| we use the mass fluxes as inputs and derive the ES effect [TTJ [32] ; and 
in section |4.3| we use the mass densities as inputs to derive asymmetric, ^-dependent 
step-edge permeability rates. In section [s] we apply perturbation theory to find P{9) 
by using primarily the mass fluxes as inputs: In section |5.1| we carry out the per- 
turbation analysis to first order for the edge-atom and kink densities as k ^ 0; in 
section [5T3] we derive the step stiffness as a function of 0; and in section [5T4] we discuss 
an alternative viewpoint on the stiffness. In section [6] we discuss our results, and 
outline possible limitations. The appendices provide derivations and proofs needed in 
the main text. 

2. Background. In this section we provide the necessary background for the 
derivation of the step stiffness. First, we describe the step configuration. Second, we 
revisit briefly the constituents of the BCF theory with focus on the GT formula and 
the step stiffness, /3. Our review provides the introduction of f3 from a kinetic rather 
than a thermodynamic perspective. Third, we describe in detail the nonequilibrium 
kinetic model [TJ |S] with emphasis on the mean-field constitutive laws for edge-atom 
and kink densities. Fourth, we set a perturbation framework for the derivation of 

2.1. step geometry and conventions. Following [71|S] we consider a simple 
cubic crystal (solid-on-solid model) with lattice spacing a and crystallographic direc- 
tions identified with the x, y and z axes of the Cartesian system. The analysis of this 
paper is for a step edge or island boundary to which there is flux / of atoms from 
the adjoining terraces. The flux / may vary along the edge, as well as in time, and 
it comes from both sides of the edge, but it is characterized by a typical size / which 
has units of (length ■ time)~^. In [T, 'E' the geometry was specialized to a step train 
with interstep distance 2L and deposition flux F, so that in steady state the flux to 
the step is / = LF. This global scenario is not necessary, however, since the analysis 
here is local and only requires a nonzero quasi-steady flux /. This could occur even 
with no deposition flux _F = 0; for example, in annealing. 

For algebraic convenience we adopt and extend the notation conventions of [S]. 
Speciflcally, we use the following symbols: (x, y, z) for dimensional spatial coordinates, 
t for time, D for any diffusion coefficient, p for number density per area , and ^ for 
number density per length; and deflne the corresponding nondimensional quantities 
5,y,5j,-D,p,|by 



(2.1) 


(i,y, S) 


= {x/a,y/a,z/a) 


(2.2) 


t 


= {af)t, 


(2.3) 


D 


= D/ia'f) , 


(2.4) 


P 


= a^p , 


(2.5) 







Now drop the tildes, so that x,y, z^t, D, p, £^ are dimensionless. This choice amounts 
to measuring all distances in units of a and all times in units of (a/)^^. Equivalently, 
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Fig. 2.1. The macroscopic (left) and microscopic (right) views of a step edge in the (high- 
symmetry) xy-plane of a crystal. In the macroscopic view, the step edge orientation relative to the 
X axis is indicated by the angle 6. The + (—) sign indicates an upper (lower) terrace. The surface 
height decreases to the right. The microscopic view shows adatoms (p), edge-atoms (<j>), left-facing 
kinks (ki) and right-facing kinks (kr); 0+ (^-) region of the upper (lower) terrace. 



(2.1)-(2.5l correspond to setting a = 1 and / = 1. For our analysis, the single most 



important dimensionless parameter is the Peclet number P from ( 1.1 1, which is equal 
to 21?^^ after nondimensionalization; i.e. 

(2.6) De = 2P-^ 

Next, we describe the coordinates of the step geometry in more detail. We con- 
sider step boundaries that stem from perturbing a straight step edge coinciding with 
a fixed axis (e.g., the x-axis). All steps are parallel to the high-symmetry ("basal"), 
xy-plane of the crystal. The projection of each edge on the basal plane is represented 
macroscopically by a smooth curve with a local tangent that forms the (signed) an- 
gle 9 with the a;-axis, where —da < < Without loss of generality we take 
< 6* < ^0 E^nd assume that 6'o < 7r/4 in our analysis. We take the upper terrace to 
be to the left of an edge so that all steps move to the right during the growth process. 



^The definition of d here is the same as that in [7], but different from the one in [S] where 9 is 
the angle formed by the local tangent and the y axis. 
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So, the projection of each step edge is represented by 

(2.7) y = Y{x,t), 

where Y(x,t) is a sufficiently differentiable function of {x,t). 

It follows that the unit normal and tangential vectors to the step boundary are |8] 

(2.8) {sme,-cose) ^ {ys,-Xs) , f = (cos 6*, sin 6*) = (x^, y^) , 

where s is the arc length and lowercase subscripts denote partial differentiation (e.g., 
Xs ■= dgx) unless it is noted or implied otherwise. The step edge curvature is 

(2.9) K = -6»3 . 

There is one more geometric relation that deserves attention. By denoting the 
densities of left- and right-facing kinks ki and kr, respectively, we have [7] 

(2.10) kr - ki = -ta,ne ; 



see section |2.3| for further discussion. This geometric relation poses a constraint on 
the total kink density, k (fc > 0). By 

(2.11) k-.^kr + h 



and (2.10), k must satisfy 



(2.12) fc>|tan( 



The formulation of a nonequilibrium kinetic step edge model (section 2.3 1 requires 
the use of several coordinate systems for an island boundary; these are described in ap- 
pendix |X] In the following analysis it becomes advantageous to use 9 as the main local 
coordinate. Its importance as a dynamic variable along a step edge is implied by the 



steady-state limit k — > | tan0| as /c ^ and P — > 0; see (4.1 1. Some useful identities 
that enable transformations to the {9,t) variables are provided in appendix [A| 

2.2. BCF model. In the standard BCF theory [6| the projection of step edges 
on the basal plane are smooth curves that move by the attachment and detachment 
of atoms due to mass conservation. The BCF model comprises the following neax- 
equilibrium evolution laws, (i) The adatom density solves the diffusion equation on 
terraces, (it) The adatom flux and density satisfy (kinetic) boundary conditions for 
atom attachment-detachment at step edges, (iii) The step velocity equals the sum 
of the adatom fluxes normal to the edge. In this setting, the GT formula links the 
normal mass flux to the step edge curvature. 

We next describe the equations of motion in the BCF model for comparisons with 
the kinetic model of section [2?3l The density, p, of adatoms on each terrace solves 



(2.13) dtp- Dt^P^F , 

where Dt is the terrace diffusion coefficient and A denotes the Laplacian in {x,y). 



As an extension of the BCF model, the boundary conditions for (2.131 are now 
formulated by linear kinetics with inclusion of both atom attachment-detachment and 
step permeability EHMHH]: 



(2.14) U=Da^{p±-P^)±DUp+-P-) 
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cf. (1.2). Here, f± is the adatom flux normal to an edge from the upper (+) or lower 
(— ) terrace, i.e.. 



(2.15) 



T f± := vp± + Dxn ■ (Vp)± 



p± is the terrace adatom density restricted to the step edge. Da is the attachment- 
detachment rate coefficient and is the permeability rate coefficient. These rates 
can account for different up- and down-step energy barriers , e.g . the ES effect in the 
case of Da^ [TT1[3^. The reference density is given by (1.3l where is replaced 
by for up- and down-step edge asymmetry. Evidently, (2.14) forms an extension 
of formula ( [1.2^ but still corresponds to near-equilibrium kinetics; it will be modified 
in section 12.31 



Equations (2.13) and (2.14) provide the fluxes f± as functions of the step edge 



position and curvature. The step velocity, v, is then determined by mass conservation. 



(2.16) 



v = f+ + f- 



In this formulation, step-edge diffusion and kink motion are neglected. In the next 
section, the BCF model is enriched with kinetic boundary conditions that account for 
the motion of edge-atoms and kinks. 

2.3. Atomistic, nonequilibrium kinetic model. In this section we revisit 
the kinetic model by Caflisch et al. [3 |8j , which is an extension of the BCF model 
(section 2.2) to nonequilibrium processes. We apply this kinetic model [7, to step 



edges of arbitrary orientation; and further revise it to account for a step edge dif- 
fusion coefficient defined along the (fixed) crystallographic x-axis. This last feature, 
although not important for our present purpose of calculating the step stiffness, ren- 
ders the model consistent with recent studies of the edge-atom migration along a step 
edge [25]. The following processes arc included, (i) Adatom diffusion on terraces. 



which is described by (2.13) of the BCF theory, and edge-atom diffusion along step 
edges, (ii) Convection of kinks on step edges with sinks and sources to account for 
conversion of terrace adatoms and edge-atoms to kinks, (iii) Constitutive laws that 
relate mass fluxes, sources for kinks and the step velocity with densities via a mean- 



field theory, and modify the BCF laws (2.14) and (2.16). In this model, kink densities 



are assumed sufficiently small, enabling the neglect of higher-order terms within the 
mean-field approach. Recently, extensions of this theory were developed [lllSl[I3]j in- 
cluding higher kink densities by Balykov and Voigt [21 [3]. Next, we state the requisite 
equations of motion in addition to (2.13) for adatom terrace diffusion. 



2.3.1. Equations of motion along step edges. An assumption inherent to 
the present model is the different kinetics of kinks and edge-atoms. Each of these 
species is of course not conserved separately, since edge-atoms can generate kinks, but 
can be described by a distinct density: (t>{x, t) for edge-atoms and k(x, t) for kinks. In 
addition, their motion is different: the edge-atom flux follows from gradients of the 
density 0; while the kink flux stems from a velocity field, w. 

We proceed to describe the equations of motion. The edge-atom number density, 
(p{x, t), solves 



(2.17) 



f+ + f- 
cos 9 



-fo 



where 15b is the step edge diffusivity defined along the high-symmetry (x-) axis and 



fo represents the loss of edge-atoms to kinks; see (2.21) and (2.27) below. For later 
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algebraic convenience, it is advantageous to transform (2.171 to {9,t) variables. By 
the formulas (A. 7) and (A. 12) of appendix [A| (2.171 is thus recast to 



(2.18) 



dt\9(j) + K{ve + vta,n9)dg(f) - Ds-^de-^de 

cos cos 



f+ + f- 



fo 



We turn our attention to kinks. The total kink density, k{x, t), of (2.11 1 solves 
(2.19) dtk + d,[w{kr - ki)] = 2{g - h) , 

where w{kr — ki) = —wtanO is the flux of kinks with respect to the a;-axis, g is the 
net gain in kink pairs due to nucleation and breakup, and h is the net loss in kink 
pairs due to creation and annihilation [7]. The terms w, g and h are described as 
functions of densities in (2.23l-(2.25l below. In the {0,t) coordinates, (2.19) reads 



(2.20) 



dt\gk + K(vg + utan 6)dek - 



COS( 



-dgiwtane) = 2(.g - h) 



Equations (2.17) and (2.19) can be transformed to other coordinates, including 
the {s,t) variables where s is the arc length. For completeness, in appendix [B] we 
provide relations that are needed in such transformations; and in appendix O we 
describe the ensuing equations of motion in the (s, t) coordinates. 



Partial differential equations (2.17) and (2.19) are coupled with the motion of step 



edges. In the following analysis, we apply the quasi-steady approximation, neglecting 



the time derivative in (2.18) and (2.20). For definiteness, the boundary conditions in 



X can be taken to be periodic. It remains to prescribe boundary conditions for atom 
attachment-detachment, i.e., specify f± in (2.15). In the present nonequilibrium 



context, f± are no longer given by (2.14) of the BCF model, as discussed next 



2.3.2. Constitutive laws. Following [7J H] we describe mean-field constitutive 
laws for fiuxes related to a tilted step edge (at 6^0). We also provide a geometric 
relation for the step edge velocity, v, which in a certain sense replaces the BCF 
law (2.16). Because the explanations are given elsewhere [2, 7], we state the mean- 



field laws without a detailed discussion of their origin. 

By mean-field theory, the terrace adatom fiux normal to the step edge is [7] 



f± = [DtP± ~ Di 



XDTP±-DK)k- 



XDtP±4' - Dxkrki) 



(2.21) 



{DTP±krki - Db)] cos 6*, j+ = 2, j_ = 3 



where I 



J' 



and Tij are (effective) coordination numbers (positive integers) that 
count the number of possible paths in the kinetic processes, weighted by the relative 
probability of a particle to be at the corresponding position. Also, Dk is the diffusion 
coefficient for an atom from a kink, and Db i& the diffusion coefficient for an atom 



from a straight edge. By neglect of Dk and Db, (2.21) readily becomes 



(2.22) 



f± = [1 + Ij^k + nij^ (p + Tij^ krki)DTP± cos 9 — De4>cqs ( 



Omitting D^ and Db is inconsistent with detailed balance, but has little effect on 
the kinetic solutions described below. 

Similarly, the mean-field kink velocity reads [7] 



(2.23) w = IiDe(I> + Dt{12P+ + hp-) - h23.DK - hDE<t) + DT{hp+ + hp-) 
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The gain in kink pairs from nuclcation and breakup involving an edge-atom is [7] 

g = cj){miDE(t) + m2DTP+ + uisDrP-) - mi2z Dxkrki 

(2.24) ~ (j){miDE(i3 + m2DTP+ + msDrP-) ■ 

The respective loss of kink pairs by atom attachment-detachment is [7] 

h = (niDEfj) + n2DTP+ + n^DTP-)krki - ^123 Db 

(2.25) ~ {niDE(j} + n2DTP+ + n3,DTP-)krki . 

In the above, 

(2.26) pij Pi + pj, pi-jk ■= Pi+Pj+Pk; p^m, n, I . 
The constitutive laws are complemented by 

(2.27) fo^wk + 2g + h , 



which enters (2.181. The step edge velocity, v, stems from a geometric relation; see 

wk + 2g + h 



appendix p3] for details. Specifically, 

fo 



(2.28) 



1 



bn cos 9 



1 



2.4. Program for step stiffness. In this section we delineate a program for the 
calculation of the step stiffness from the model of section 2.3 The key idea is to reduce 



the nonequilibrium law (2.21 1 to the linear kinetic law (2.141 by treating the normal 
fluxes, /±, as external, free to vary, 0(1) parameters of the equations of motion along 
a step edge. In this context, the diffusion equation (2.13 1 is not invoked 

i and k 



Our method 
The solution studied 



relies on the perturbation of a solution for the densities 
here is that of the kinetic steady state, under the assumption that it can be reached. 
Accordingly, we neglect the time derivative in the zeroth-order equations of motion; 
furthermore, we neglect this derivative to the next higher order by imposing the quasi- 
steady approximation. Another case, left for future work, is that of thermodynamic 
equilibrium; see section |6] In summary, we apply the following procedure: 

(i) To extract the kinetic steady state, we set dt\e = and k = (i.e., we consider 
straight edges). This leads to a system of algebraic equations for {(j), k) = (0^°\ fc(°))[^ 
The coefficients of this system depend on 9 and f±. In principle, (</'''*'■', fc'^'^') cannot 
be found in simple closed form at this stage. 

(ii) We assume that P ^ 1, and determine relatively simple expansions for 
(0("), fc(o)) in powers of P for < 61 < 0{P^/^) and 0{P^I^) < 9 < tt/A . 

(iii) We replace (0, fc) by (0(°\fc(°)) in the constitutive law (2.211 and compare 
the result to (2.141. Here, our analysis follows up two mathematically equivalent but 
physically distinct routes, (a) By taking f± as input parameters, we derive formulas 
for the adatom reference densities, pf, and attachment-detachment rates, D^^, that 
depend on f±; cf. (2.141. Step permeability is not manifested in this setting {Dp = 0). 
(b) By considering p± as inputs, we predict attachment-detachment rates and non- 
vanishing step permeability rates. 

(iv) We consider perturbations of the kinetic steady state by taking < |k| ^ 1, 
i.e. slightly curved step edges. Accordingly, we let 



(2.29) 



^In this context, the superscript in parentheses denotes the perturbation order in k. 
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where 



and k/c*-^-* are d eviations from the kinetic steady state and depend on 



{(j)^^'> , fc(°^ ) . Expansion ( 2.29 1 is imposed on physical rather than mathematical grounds. 



Indeed, if the mean-field flux (2.211 is expected to reduce to the linear kinetic law (2.14), 
then (p must be linear in k. The equations of motion along an edge and the constitutive 
laws are linearized in K(j)^^^ and Kk^^\ 

(v) By treating f± as input external parameters, we replace 4> and k in the 



right-hand side of the constitutive law (2.221 by expa nsions (2.291. Su bsequ ently, we 
determine the stiffness I3{9; /+, /_) by comparison to (2.141 in view of (1.3). 

The choice of fluxes /± or densities p± as input parameters is a physics modeling 
question. Although the mathematical results are equivalent for the two choice, the 
physical interpretation of these results is different, as stated above. 

3. Main results. Here, we give the main formulas stemming from our analysis of 



the kinetic model described in section 2.3 A necessary condition for our perturbation 



analysis is < < 0{P) ^ 1, to be shown via a plausibility argument in section 5.1 
Derivations and other related details are provided in sections |4] and |5] 



3.1. ES effect (section 4.2). When the fluxes f± are input parameters, the 
attachment-detachment of adatoms from a terrace to an edge is asymmetric. So, the 
related diffusion coefficients Da^ , or attachment and detachment kinetic rates, which 



enter (2.14), are found to be different for an upper and lower terrace: 

■i2 



(3.1) 



Da+ 

Da' 



Dt[1 
Dt[1 



TO2I 



TO31 



/,(0) 



\n2 



tan 9) cos ( 



tan"' 



where < 9 < 7r/4 and (^2, '^■2) 7^ (^3, m^a, "■a)- For < P <C 1, we show that (3.1) 
reduce to 

(3.2) Da"^ Dt{1 + lj± tan 6*) cos 6* , 

where j+ — 2 and j_ — 3. In this description, there is no step permeability. Note 
that the results presented in this section and their derivations do not depend on the 
step edge curvature. 

3.2. Step permeability (section 4.3). By using the adatom densities p± as 



inp ut ext ernal parameters, we show that step permeability coexists with the ES effect; 
cf. (2.14). For 0(P^/^) < 6 < it /A the diffusion coefficients for permeability are 

Azr(l + L- , tan( 

(3.3) " 



= Dt 



■ cos 



1 + {A+ + A_)cos( 

The accompanying (asymmetric) attachment-detachment diffusion coefficients are 

1 + Ij^ tan6' ± A^{l2 - h) sin( 



(3.4) 



where 



(3.5) 



Da^ =Dt- 



tan6' ± A^{l2 
1 + {A+ + A_)cos9 



■ cos 9 



1 1 -f ^3 tan 9 



Q(l) 



A^ = 



1 1 + /a tan 6* 



sin 61 Q(l) 



(3.6) Q(p) =pi(l + /2tan6i)(l + ;3tan6l)+p2(l + ^3tan6')+p3(l + ;2tan6l) , 

with p (pi, p2, P3) and p — I, m, n; in (3.5), 1 = {li, I2, h)- Note that 
and Da^ here are independent of f±, as in (3.2). The corresponding results for 
< 9 < 0{P^/^) are presented in section 4.3 Again, the results presented in this 



section and their derivations do not depend on the step edge curvature. 
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3.3. Step stiffness (section 5.3). Let the adatom fluxes f± from an upper 
(+) and lower (— ) terrace towards an edge be the input, independent parameters. 
For sufliciently small angle 9, the stiffness is found to be 



(3.7) 



(3.8) 



kbT ni23 J+ + f- V 



/3 



p-2/3 



4^1 



23 



(9(p-2/3) < 6* < 0(P^/^) < 1 



ksT ' ?^123(C'o)^ + 8^123 

where P123 (p = m, n) is defined in (2.26) and[^ 

2toi23 



(3.9) 



"-123 ' 



123 



-(/+ + /-) 



1/3 



Matching the asymptotic results (3.7 1 and (3.8 1 is discussed near the end of section 

EH 

For 9 = 0(1) the formula for P becomes more complicated; we give it here for 
completeness. Generally, 



(3.10) 



where (j)^^^ and 0^^^ are expansion coefficients for and depend on f± and their 

ie coefficients are obtained 

'4, 

C^P, (j)^^'>^CfP, 



derivatives in 9; cf. (2.291. These coefficients are obtained explicitly for P ^ 1. In 
particular, for 0{P^'^)< 9 < 7r/4, 



(3.11) 



(3.12) 



. _ 1 (l + ^3tang)/+ + (l + /2 tang)/, 
^0 " 2sin^? 



(3.13) cf 



(3.14) 



Q(l) 



(4"^ + tan 9)k''°^ cos 9 + (w^o) tan 9)e tan 9 + w^°^ + 



2 sin 61 



/o 1 ir\ (0) r,7 /~«t> , 1 ^xf+ + Sa;/_ + 2C^ 

(3.15 > ^2liC^ + I2——— — ^ + ^3 1 , , , — ^, S:, = (cos6') 

1 + t2 tan e' 1 + t3 tan 9 

(3.16) Vr"^ = 2/i 



2^2 



2^3 



1 + ?2 tan 6* 1 + U tan 9 ' 



(3.17) M/'^' = -h 



h + f tan ( 
(1 + I2 tan 6*; 



Ks./+ + 2Co^)-Z3 



/3 + tan 6* 
(l + /3tan6')2 



(s,/_ + 2C*) , 



•^The superscripts in C'', C"^, C*" and elsewhere below indicate the physical origin of these 
coefficients, and should not be confused with numerical exponents or perturbation orders. 
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(3.18) 



H" = 



tan ( 



s,f+ + 2Cl 
l + h tan 9 



"3 



\ + h tan 9 



(3.19) 



tan 9 



are coordination numbers. Recall definition (3.6 1 for Q(l). It is 



where rij and Ij 

worthwhile noting that there is no asymmetry in the step stiffness /3, in contrast 
to the attachment-detachment coefficients. The reason for this difference is that $ 



depends only on the edge-atom density, as shown in (3.101 



For the alternative approach in which the adatom densities p± are specified rather 



than the fluxes /±, the analysis of the step stiffness is presented in section 5.4 The 
corresponding result (5.44) is not of the form (1.2) and (1.3), however, since the 



coefficient f5 in (5.45) is not proportional to p*. 



4. The kinetic steady state . We analyze the kinetic steady state for a straight 
step, including its dependence on the Peclet number P in section [41] and the ES effect 
and step permeabiity in sections |4.2| and |4.3[ 



4.1. Kinetic steady state and its dependence on P. In this section, we 
simplify the equations of motion for edge-atom and kink densities by imposing the 
kinetic steady state {dt = 0) for straight steps (k = 0). We find closed-form solutions 
for small Peclet number, P ^ 1, in two distinct ranges of 9. For 9^ = 0(Fi/3) <9 < 
7r/4, we show that <\) — <\)^^^ is given by (3.11 1, and k = fc(°^ is given by (3.19), or more 
precisely by 



(4.1) 



fc(°) - tan 61 + C^P, 



where 



(4.2) 



_ 2Cl 2ClQ{ra) cosg m^jl + ^3 tan6>)/+ H- m-i{\ -f h tan6>)/_ 
° ^ tan 61 2C;fQ(n)cos6l-hn2(l-H?3tan6')/+-t-n3(l + /2tan6i)/_ 



(5(p) and Cq are defined by (3.6) and (3.12). Furthermore 
(4.3) ~ Ct P'/' , - Cl 



< 6* < 6*^ = 0{P^I^) 



where is defined by (3.9), 
(4.4) 



123 



and pi23 {p — I, m, n) is given in (2.26); cf. equations (4.27) and (4.28) in 8 . In 
effect, we determine mesoscopic kinetic rates, including the attachment-detachment 



and permeability coefficients in (3.1l-(3.5) 



We proceed to describing the derivations. By dt\g 
and (2.28), we have f+ + f- ~ fncos9. q = h and v 
terms of De4> ~ 2P 



and K = in (2.17), (2.191 



/pcosff. Elim inate Dtp in 
using (2.22). Thus, we readily obtain (3.14) for u^^^ := v, 
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along with the following system of coupled algebraic equations: 

[mi0(°) - (fc(o)' - tan2 6)] [l + hk^'^^ + m20(o) + =f - tan^ 9)] 

X [1 + ?3fc(o) + m30(°) + ^|i(fc(o)' - tan2 0)]2P- + K^/.^") - ^^(fc(o)' - tan^ 6)] 

X [1 + ?3fc(") + m3,^(°) + - tan2 e)] {s^f+ + 2P- 

+ [m30(o) - ^|i(/c(0)^ - tan2 0)] [l + l^k^^) + ma^f") + f (fc(o)^ - tan^ 9)] 

(4.5) x(s,/_ + 2P-V^°^) = , 



- (fc(o) 



2 



"2 I 



aifc(°)+3mi(/.("))[l + Z2/fc(°^ 
X [1 + ;3fc(") + m3 
x[l + ?3fc(°) + m3 
X [1 + ?2fc^°^ + m20(°) + (fc(o)^ - tan2 0)] {s^f_ + 2P- Vo) 
= (/+ + /-) [1 + ^2^0 + m20o + f (fc(°)' - tan2 9)] 



-tan2 0)] 
tan2 0)]2P-V^°^ + (?2A:<°^ 
tan2 0)](s,/+ + 2P-V(°)) 



3m20(°)) 



(4.6) x[i + ;3fc^*'^ + 



TO3 



(fc(o) ~tan2 0)] 



Once these equations are solved, the flux variables w —: g —: g^^^ and h =: 
are determined in terms of f± by the constitutive laws ( 2.23 )-( 2.25 1 . The substitution 
of (j) and k into (2.221 provides a relation between f± and p±. 

Next, we simplify and explicitly solve (4.51 and (4.6 1 by enforcing P ^ 1. The 
ensuing scaling of < 



and fc*^*^' with P depends on the range of 9. We distinguish the 
cases Oc{P) < 9 < it / A and < 6* < 9c{P), where 9c is estimated below; we expect 
that 6*^ ^ as P 0. 

(i) 9 = 0{1). By seeking solutions th at a re regular at P = 0, we observe that if 
P = then (^(°),fc(°)) = (0,tan6') solves (ITsl and Thus, the expansions 



(4.7) -C^P , fc'"^ - tan + P , Cf" = 0(1) 

form a reasonable starting point. These expansions yield the simplified system 

2C^(2miC^ - niC^'tan6l)(l + Z2tan6')(l + I3 ta.n9) 
+{2m2C^ - naC^ tan 6')(1 + ^3 tan 6l)(s^/+ + 2C^) 

(4.8) +(2m3Co* - n3Co^tan0)(l + ?2tan0)(s,/_ + 2Ct) = , 



2C;f /itan6'(l 



(4.9) =s, 



+?2 tan 61(1 

:(/+ 



h tan6')(l + h tan( 
^3tan6')(s2,/+ 



2Co*) + /3tan6'(l 
/_)(l + ;2tan6l)(l + /3tan6l) . 



-h tan0)(s,/_ + 2Co^) 



The solution of this system leads to (3.121 and (4.2| 



We now sketch an order-of-magnitude estimate for 9c, the lower bound for 9 in 
the present range of interest. By (4.2 1, Cq = 0{\/9'^) for 6'c < 6* < 1. Hence, 
expansion (4.7) for the kink density fc^-* breaks down when its leading-order term, 
tan 6*, is compar able t o th e correction term, C^P: 9c — 0{P/9^) by which 9c = 
0(pi/3). Thus, (|4j|-(|49| hold if 0(pi/3) < 61 < 7r/4. A more accurate estimate of 



ANISOTROPIC STIFFNESS FROM EPITAXIAL KINETIC MODEL 



15 



the lower bound requires the detailed solution of and MM for 9 = 0{P^^^), and 
will not be pursued here. 



(ii) < 6 < 0(P^/^). For all practical purposes we set = in (4.5 1 and (4.6 1 
We enforce the expansions 



(4.10) 



Ct''" - 0(1) as P ^ 



and find the exponents v and a by reductio ad ahsurdum. The only values consistent 



with (4.5) and (4.6) readily turn out to be 

(4.11) y = 2/i, 



1/3. 



These values are in agreement with the analysis in [H]. By dominant-balance argu- 
ments, the coefficients Cq and Cq satisfy 



(4.12) 



2^123 Cq Cq — f+ + f- 



by which we readily obtain (3.9) and (4.4). Note that the zeroth-order kink velocity 
becomes 



(4.13) 



,(0) 



2h23P-''^Ct = 0{P-'/') 



(iii) Consistency of asymptotics for 9 — 0{P^^^). As a check on the consistency 
of our asymptotics and the estimate of Oc, we study the limits of (4.3) and (4.7) in 
the transition region, as 9 ^ 0(P^/^). It is expected that the two sets of formulas for 
(I) and k, in 0(P^/^) < 9 < n/A and < 9 < 0(P^/^), should furnish the same order 
of magnitudes. 

Indeed, by letting g -> 0{P^^^) < 1 in (jZ?) for 
in agreement with (4.3) for — 2/3 



we find = 0{P/9) ^ OfP^/^), 

Similarly, setting 6* = 0{P^/^) in (Ej) for k 



5.3 



yields k = 0{9) 0{P^^^), which is consistent with (4.3) for a = 1/3. In section 
we show that such a "matching" is not always achieved for the first-order corrections 
0*^^^ and k'^^\ since the corresponding asymptotic formulas involve derivatives in 9. A 
sufficient condition on the ^-behavior of the fluxes f± is sought in the latter case. 



In the following, we use the kinetic steady state in the mean-field law (2.211 



to derive mesoscopic kinetic rates as functions of 9 by comparison to the BCF-type 
equation (2.14). We adopt two approaches. In the first approach, f± are used as 



external, input parameters; the effective kinetic coefficients are thus allowed to depend 
on f± . In the second approach, the densities p± are the primary variables instead. 

4.2. Flux-driven kinetics approach: ES effect. In this subsection we treat 



the fluxes f± as given, input parameters. Accordingly, we derive (3.1) and (3.2), i.e. 
the attachment-detachment rates -D^^ for adatoms. 
reference densities pf entering (1.3) are 



In addition, we show that the 



(4.14) 



2Ct 



DTil + lj± tan( 



j+ - 2, j_ - 3 , 0(pi/3) <9<7r/A 



(4.15) 



2c; 



D 



f^p-i/3 o<e< o{p'/') , 



where Cg and Cq are defined by (3.12) and (4.4|. 
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(4.16) f._ 



The relevant derivations follow. The substitution of (4.7 1 into (2.221 yields 

l + L-^fc(°^+TO. 



"|^(/fc(0)'-tan2 



DtP± cos0-2P~^ 



COS! 



Here, we view the linear-in-p± term of (4.161 as the only physical contribution of the 
ad atom densities to the mass flux towards an edge. Consequently, by comparison 
to (2.141, the coefficient of this term must be identified with Da^- Thus, we extract 
formulas (3.1 ). In addition, we obtain = 0; so, step permeability is not manifested 
in this context. The reference density of ( |1.3[ ) is 

(4.17) "± 2P-V(°)cosr 



Da' 



which is in principle different for an up- and down-step edge. 

Note that Da^ and pf^ depend on f± within this approach. Further, the ratio 
of Da^ and Da~ depends on the values of Ij, nij and nj. For suitable coordination 
numbers, it is possible to have Da~^ > Da~ , i.e. a negative (vs. positive) ES effect [TTJ 
132] , which can lead to instabilities in the step motion. Next, we derive simplified, 
explicit formulas for Da^ and pf when P <C 1. 



(i) 0(pi/^) <e < 7r/4. By substitution of (|4j| with ( |3.12| and (|4^ into ( pU l 
we have 



(4.18) 



Da+ 
Da' 



Dt [1 + h tan e+{l2+ \n2 tan 9) C^P + msC^P] cos 9 , 
Dt [1 + ^3 tan 61 + (^3 + ^713 tan 9)C^P + mgC^P] cos 9 , 
which reduce to (3.2) as P ^ 0. 



In the same vein, by (4.17) the reference densities are 

2Cn* 



(4.19) 



* DT[l + l2tB.n9+ {h + ^n2ta.n 9) C^P + m2C^P] ' 

< 

* Dt [1 + ^3 tan 61 + (?3 + in3 tan 61) C^-p + maCo^P] ' 

which readily yield ( 4.14[ ). 

(ii) < 9 < 0(P^/^). In this case, we resort to (4.3). Equation (3.1) for the 
kinetic rates furnishes 



Da 
Da 



(4.20) 



Dt{1 + hC^ P'/' + [m2Ct + \n2{Clf]P^'^]cos9 
Dril + hC^ P'^'')cos9 , 

Dt{1 + hC^ pi/3 + [m,C^ + in2(Co^)2]p2/3} COS0 
i:»T(l + Z3C'o ^'^^^)cos6' . 



To leading order in P, these formulas connect smoothly with (4.18) and, thus, jus- 
tify (3.2) for < < 7r/4. Furthermore, pf are given by 



2C^ p-i/3 



2c; 



(4.21) 



Dt{1 + hC^ P^/'') 
2Ct P-^/^ 

~W + hdfp^) 



Dt 
2Cf: 



Di 



p-^'^i-hc^p^''') 

p-l/3(l_;3C^o*pl/3) 
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which reduce to (4.151. Notably, pf depe nd o n the fluxes, f±, through Cq. 

A few remarks are in order. First, by (3.2) the ES effect is present for 0(P^/^) < 
9 < tt/A only if I2 7^ h- Accordingly, our formalism provides explicitly an analyti- 
cal relation between the number of transition paths for atomistic processes and the 



mesoscopic kinetic rates. Second, formulas (4.201 show that for P <C 1 the nonzero 
ES barrier is a corrective, 0{P^^^) effect for sufficiently small 9, even when I2 7^ ^3. 

4.3. Density- driven approach: Step permeability and ES effect. In this 
subsection we show that the treatment of the densities p± as independent, external 



parameters in the kinetic law (2.221 leads to coexistence of the ES effect and step 
permeability. In particular, the permeability and attachment-detachment rates are 



provided by (3.3)-(3.5) 



To derive (3.3) and (3.4), we solve (4.16) for f±, which are viewed as dependent 
variables, taking into account that tf>^^'> and fc'^'^^ depend on f±. To simplify the algebra 
while keeping the essential physics intact, we restrict attention to 0{P^^"^) < 9 < 7r/4. 

First, in view of (4.7 1 we further simplify relation (2.22). By 



(4.22) 



where A± are defined by (3.5), the adatom fluxes at the step edge reduce to 



(4.23) /± (1 + lj± tan 9)DtP± cos 9- (A+/+ + A_/_) cos( 



P< 1 



Second, we invert (4.23) to obtain f± in terms of p±. Equation (4.23) reads 



(1 + A+cos6')/+ +cos6'yl_/_ = (1 + /2 tan 6')i:'TP+ cos 6* , 
cos9A+f+ + {l+cos9A^)f^ = (1 + tan 6')DtP- cos 6* . 



(4.24) 

The inversion of this system yields 
■(1 + /2 tan 6*) ( 1 -I- cos 6* ^_) 



/-f 
f- 



1 + {A+ + A_)cos( 
(1 + I2 ta.n9)A+ cos 6 
1 + {A+ + A^)cos9 



Dtp 



DtP+ 



A_(l + Zatan^) cos 6* 
^ " 1 + {A+ + A^)cos9 
(1 + A+ cos0)(l + ?3tan( 



1 + {A+ + A_)cos( 



Dtp- 
-Dtp 



These relations have the form of the kinetic law (2.14); by comparison, the rates D: 



are given by (3.3 1, while 
(4.25) 



4 







This value is expected since the system is homogeneous in this setting, i.e. f± =0 only 
if p± — 0. The reference density pQ becomes nonzero (but small in an appropriate 
sense) if we allow in the formulation nonzero values for Db and Dk, i.e. nonzero 
diffusion coefficients for an atom to hop from a kink and a straight edge. The study 



of these effects lies beyond our present scope. Equation (4.25 1 challenges the definition 



of the step stiffness; see section 5.4 



Equations (4.24) also predict an ES effect. Indeed, by recourse to (2.14), the 
related attachment-detachment rates are 

, (1 + L.tan9)(l + A^cos9) . 
(4.26) P'a = ^ ^"^^ — , , , ,^ — -Dt cos 9-D^ 



1 + {A+ + A_)cos( 



which readily yields (3.4) by use of (3.3). 
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The behavior of the fluxes 
< 9 < 0(pi/3). Indeed, by (|4.3|)-(|4.4| the density 
function of /+ + 



as fun ctions of p± is dramatically different for 
is a nonlinear algebraic 
/_ in this case. Thus, the mean-field constitutive equations in 
principle cannot reduce to kinetic laws that are linear in p± . This approach does not 
lead to standard BCF-type conditions at a high-symmetry step edge orientation. The 
implications of this behavior warrant further studies. 

In the following analysis for the stiffness we emphasize the flux-driven approach. 

5. Perturbation theory and step stiffness. In this section we consider slightly 
curved step edges, and apply perturbation theory to find approximately the edge-atom 
and kink densities, </) and k . from the kinetic model of section |2.3| On the basis of 
the linear kinetic law (2.141 along with (1.3 1 for po, we calculate the step stiffness, /3, 
as a function of the orientation angle, 9; see formulas (3.7l-(3.19). The underlying 
perturbation scheme for the densities is outlined in appendix [E| 

The starting point is expansion (2.291, which we assume to be valid for 
9 < tt/4: and view as a Taylor series, 
the kinetic steady state of section 



4.1 



< 

The functions 0'-°' and fc'"' correspond to 

The first-order coefficients (f)''^^ and k'^^^ are 

locally bounded and are evaluated bel ow. O nly the coefficient (j)^^^ is needed for the 
calculation of the step stiffness, f3, by (2.22); for completeness, we also derive k^^\ 

The relation of (3 to 0^°) and (j)^^' is provided by the following argument. By 
substitution of (2.291 into ( |2.22[ ) and treatment of f± as given external parameters 
(in the spirit of section 4.2 1, we obtain 



(5.1) 



= [l + ;,±fc(")+m,^0(") + ^(fc(o)'-tan20)]i?7.p±-i?i;0(o)-Ki?B^(i) 



cos 



where j+ — 2 and ~ 3. By comparison of (|5.1|) 
2P-1 = De) 



to (1.3) and (2.14), we have (using 



(5.2) 



Da^P: 



A. 

knT 



2P- 



-1 



by which we assert (3.10) in view of (4.17). Our task is to calculate c/)'"'^' in terms of 



and P when P ^ 1. 



5.1. Linear pert urbat i ons. In this subsection we derive formula (3.11 1 for 0^^) 
along with (|3.13| and (|3.16|)~(|318| when 0(pi/3) < < 7^/4. in addition, we show 



that in this regime 



(5.3) 



^(1) {v^°^ + v^°h'An9){cos9)~^ + {w^"h&n9)e 



For < 6* < 0(pi/3), and fc^^^ are 
(5.4) 0(1) ~ 4^123 



2^?*= cos 61 



= 0(1) , 



(5.5) 



-AP-^/^h 



23 



^0 



ni23{C^r + 8^123^0^ 



Recall that Cq and Cq are defined by (3.9) and (4.4). Furthermore, we demonstrate 
that \k\ should be bounded by P for the perturbation theory to hold; see (5.40). 
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We proceed to carry out the derivations. Following appendix |Ej we formulate a 
2x2 system of linear perturbations for and k. First, we linearize the algebraic, con- 
stitutive laws (2.23 H p. 25 1. Expansions ( |2.29 1 induce the approximations w{(j),k) ~ 
+ ~ g{(l)^°\k^°'>) + Kg(^^ and /i(<^, A:) ~ h{<l)^°\k^°^) + Kh^^\ 

where 



(5.6) 



In addition, = f+ + f- and 5(0) := g((/)(o), = fa((/)( "), fcW ) h^°l Second, 
we replace the above expansions in the equations of motion (2.181 and (2.201 and the 



constitutive law (2.271. Hence, we find the system 



(«;^fc(°) + 25^ + + {wkk^°^ + + 2gk + hk)k^^^ = -(^"^ + v^°^ tan( 

(5.7) 2{g^ - /i^)0(i) + 2(5fc - hk)k<^^^ = (1;^°^ + tan 0)4°^ + tan0)e , 

where w'^"^ := w((/)'^*'^, fc'^*'^) and Sx ~ 1/ cos0. This system has solution 



(5.8) 
where 
(5.9) 
(5.10) 



2?0 



-(w^"^ + wW tan 61)0, 



+ {2g + h)^ Wkk(°^ + + {2g + h)^ 



2{g~h)^ 



2{g-h)k 



(0) 



{vf^ + tan6')4°^ + S:,(u'('') tan6')e 2{g - /i)fe 



Wfcfc(o) + 2(g + /i)fc 



a(0) 



(5.11) 2?'= 



w^fc(o) + (25 + h)^ -{vf^ + tan e)(j)f^ 

2{g - h)^ (1;^°^ + tan^)^;^") + s^(t«(°) tane)^ 



Note that 0'^-'^) and k^^^ depend on the ^-derivatives of the zeroth-order (kinetic steady- 
state) solutions. 

By ( |2.23 H( 2.25[ ), we calculate the 0- and fc-derivatives of w, g and h: 

( 2P^^ 

1 1 + ?j,fc(o) + mj,0(o) + ^(fc(o)' - tan2 0) 

(cos0)-i/, + 2P-VW \ ■ 9 ■ , 

(5.12) — TO," ^ — s > , 7 1 = 2, 7_ = 3 , 

[1 + A:(0) + m, (^(0) + ^(fc(")' - tan^ e)Y I 



(5.13) = - E + 



(cos6')-i/g + 2P-V^°^ 



[1 + ;j,/c(o) + mj,0(o) "f (fc(o)^ - tan2 6')]2 ' 



ff0 



f>=^ — ^1 



(cos0)-i/9+4P- 



+ l^m + mj^0(o) -f ^(fc(")^ - tan2 



(5.14) -mj,(/) 



(0). 



(cOS0)-l/g + 2P- 



[1 ?j fc(o) TO^. 0(0) + !^(A:(0)2 _ tan2 6l)]2 
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(5.15) 



^ [1 + Ij m + mj 0(0) + ^(fc(o)^ - tan2 e)f 



q=+ 



2P-in 



1 + I J m + TOj 0(0) + ^(/c(o) - tan2 



(5.16) -mj^Uj^ 



(cos6l)-i/g + 2P-V^°^ 
[1 + ?j,fc(") + mj,0(o) + ^(fc(o)^ - tan2 6l)]2 



fc(0) 



E 



(5.17) -l(fc(")'-tan2 0) ^ 



(cos6i)-Vg + 2P-V^°^ 
1 + l^m + mj,0(") + ^(fc(o)^ - tan^ 9) 

(Z,, + ^fc(o)) [(cos0)-i/, + 2P-V(°)] 



[1 + fc(o) + TUj 0(0) + ^(fc(o)^ - tan2 6l)]2 



Equations (5.8h(5.17| are simplified under the condition P <^\, which we apply 



next. We distinguish two ranges for the angle 6. 

(i) 0(pi/3) <_e< 11 /A. We proceed to show ( |3.1l| and ( |3.13[ ) for 0(i). By 
using (4.7l with (3.121 and (4.2), we replace 0(o) and by their expansions in P. 



Thus, the derivatives of w, g and h are simplified to 

2h , 2/3 



(5.18) u;^~F-M2?i + - 



^2 tan 9 1 + ^3 tan 9 



P-^W^ = O(P^i) , 



(5.19) u.,^-Y. h,^^. + l^tan0) (^°:'^"i+,',fi = 0(1) , 



</=+,■ 



(1 + ;^, tan 61) 2 



(5.20) 



, (cos^)-V, + 4Co^ 



= 0(1) 



9=+,- 



(5.21) gk^-PCt V m,,(Z,, + >tan( 



{cos9)-^f^ + 2Ct 
(1 + /^, tan 6*) 2 



= 0(P) 



"2 



"^3 



(5.22) /irf, ~ C^' tansfni + , , , , . 

^ ' ^ ° \ 1 + ^2 tan 61 1 + /a tan 61 



0(1) , 



tan 6* 



(5.23) h,^'-^i2mCt+ 



icos9y% + 2Cl 
1 + Ij^ tan 



= 0(1) 



It follows that the determinants of (|5.9|-(|5.11| are 
(5.24) 



V 2P-^ hkW'l' tun 9 
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(5.25) 2?"^ ~ -{wk ta.n0 + w^°'> + /ife)[(4°^ + v^°'> tan 6')fc^°^ + (cos e)-\w^°^ tan 0)g] , 



(5.26) 



pfc ^ P-1W^**^[(t;(°) +t,W tan0)(cos( 
cos & 



+ (w(°) tan 



Hence , in view of (5.8l, the coefficient (/>(^) is given by (3.11 1 with (3.13) and ( 3.14 )- 
(3.19|) under the replacements H'' := hk, W'l' := Pw^ and VF'' := Wk- By (|5.8|), the 



corresponding coefficient fc^^^ is given by (5.3) 



(ii) <9 < 0{P^/^). We now calculate the first-order corrections 0^^) and fc'^ 



by (^5.8|-(5.17) with recourse to formula (4.3) with (3.9) and (4.4| 



We start with (5.8 1. The requisite derivatives of w, g and h in the present case 
(where practically 6 = 0) reduce to 



(5.27) 
(5.28) 
(5.29) 
(5.30) 
(5.31) 
(5.32) 



gk 



W4, ^ 2Zi23P-^ = O(P-I) , 
w,^-2P-^l\ll + ll)Ct = 0{P''''), 

^ 4rm23P-'^' = 0{P-'/') , 
^ -2pi/3(m2/2 + m3h){Ctf - 0(pi/3) ^ 
~ hni23P-'/HC'or = 0{P-'/') , 



Note that w^^^ is given by (4.13) 



It follows that the determinants V, and V'' of (5.9l-(5.11 1 become 
(5.33) 



(5.34) 



V ^ -P-'/''C^h23[ni23{C^)' + 8mi23C'o^) = 0{p-''/') 



(5.35) 



Since dg{C^) is finite at 6' = 0, (5.4) and (5.5) ensue directly via (5.8). 

(iii) Transition region, e = 0(pi/3). Next, we study the limits of the and 
fc(^' found above when 9 enters the transition region, 9 0{P^^^). 



„ _ . _ „ ^ ^ „ " \ / 

First , we co nsider (f)'''^^ in the range 6 > 0(P^/^) and take 6* <C 1. By (3.12) 
and ([3!6t-( |3.19[ ), we find H'' = 0(1), W'' = 0(1/6'), W'>' = 0(1), and 



(5.36) w 



{w^°^0)e = if+ + f-)e\e=o + 0{9) 



(5.37) {v^°^ + tan0)fc^°^ + {cos 9)-\w^^^ tane)^ = 2(/+ + /-)e|e=o + 0(0) 
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Hence, assuming (/+ + /-)e 7^ at 9 = 0, we have 

(5.38) cl)^^'> =O{P/0^)-O{{f+ + f-)e) 0{P^/^)<e^l 



which becomes 0(P^/^(/++/_)6/) as ^ ^ 0(P^/^). On the other hand, by (5.4) we get 
0(1) = 0(1) when 0<9 < 0{P^/^). This behavior is not in agreement with ( |5.38[ ) un- 
less if++f-)e = 0{P-^'^), i.e. the fluxes vary over angles 0(pi/3)^ ^ f±{P^^'^0) 
for 6 = 0(P^/^). This behavior of f± is not compelling, since it is generally expected 
that the agreement in orders o f ma gnitude is spoiled by the ^-differentiation. 

We next consider k^^\ By O we find fc^^) = 0((/++/-)e) for 0{P^/'^) < 6I < 1. 
On the other hand, by fc^= 0{P~^/^) for < 6* < o\p^/'^). The two orders 

of magnitude agree if (/+ + f^)e = 0{P-^'^) as above. 

5.2. Condition on k, and P. Thus far, we have not provided any condition for 
the validity of our perturbation analysis. Such a condition would impose a constraint 
on K, and P. In principle, k is a dynamic variable. For appropriate initial data, 
the step edges are assumed to evolve to the kinetic steady state with k = 0. Small 
deviations from this state can be treated within our perturbation framework if 

(5.39) |Kfc(i)| <A:(o) . 

By revisiting the formulas of sections 4.1 andjs] for and k^J\ we can give an 
order-of-magnitude estimat e of an upper boun d for n. By comparison of the 0{P) 



correction term for fc(°) in (4.11 to fc^^^ in (5.3 1, where 9 = 0(1), we obtain 



(5.40) 



5.3. Step stiffness. Once 



k\ < 0{P) . 

1) have been derived, the step stiffness 



and 



follows. We invoke the formulation of section 5.1 on the basis of formula (3.10) by 



using the fluxes f± as input external parameters. In particular, we show the limiting 

for small 



behaviors (3.7) and (3 



In correspondence to section 5.1 



we use two 



distinct regimes. 

(i)O(P^) < 9 < 7r/4. By (5.2) and the analysis in section 
by ( pJ0| -( |3l9| . Speciflcally, 



5.1 



(3 is given 



(5.41) 



13 Of 



keT Ct ' 



which is an 0(1) quantity in P when 9 = 0(1). In order to compare this result to a 
recent equihbrium-based calculation for the stiffness [37] i we take 0(P^/^) < <C 1. 
Then, by pl2| . 



(5.42) 



f «123 



In addition, if (/+ + /_)e 7^ as 6* ^ 0+, by (3.13) and (5.37) we flnd 



(5.43) 



CI 



4, (/+ + f-)9\e=o 



ni23 9^ 



O 



Thus, (|3.7[) follows from (5.41 ). By contrast, if (/+ + /-)6/ vanishes in the limit 9^0 



then, by (5.37), l3/{kBT) = 0(1). 
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(ii) < 61 < 0(pi/3). In view of ( |3.10| with (jTsf and ([H^, we readily obtain 
formula (3.8 1 for 13. 

(iii) e 0{P^/^). Formula is consistent with the 0(1/6') behavior of 
for 0(P^/^ ) < 61 < 1 provided that (/+ + = 0(P-i/3). Indeed, from KIl ) 
via (5.371 we have P/lksT) — 0{{f+ + f^)g/6), which properly reduces to (3.8l. 
Again, this "matching" is not compelling since ^-derivatives are involved. 

5.4. Alternative view. We consider 6 = 0{1) and focus briefly on the implica- 
tions for the stiffness of treating the adatom densities p± as input parameters. This 
approach is mathematically equivalent to that of section |5.3[ only the physical defini- 
tions are altered in recognition of p± as the driving parameters. This viewpoint was 
partly followed in section 4.3 for straight step edges (k = 0). 



We show that the adatom fluxes have the form 



(5.44) f±^DA^p±±D^ip+-p^)~&{e;P+,P~)-^ ■ 

The coefficients Dp{9) and Da^{0) are defined by (3.3) and (3.4|; and 



(5.45) 



2Cf 
{A+ + ^_)cos( 



cose o{p^''^) <e <TTi 



where Cf and A± are defined by (3.131 and (3.51. Furthermore, the /^-dependent 
Cf is now evaluated at f± 
Notably, 



Da P± ± -Dp (jO-|_ — (0_); thus, 6 becomes p-dependent. 



(5.46) 



6 = 0(1/61) 0{P^^^) <e <t:i 



As noted in section |3.3[ these results do not have the usual form since 6 is not pro- 
portional to /o*. 



We derive ( 5.44 )-( 5.46 1 di rectl y from (5.1 ) by treating the term n 0^^) as a pertur- 
bation. For Kcj)'-'-^ = (section 4.3l, (2.141 for f± is recovered with po = 0; see (4.251 



For K ^ 0, (5.1 1 reads 



(5.47) 



/± - (1 -f ta.TLe)DTP± cose - {A+f+ + A_/_) cos6l - Cf kcos( 



By viewing Cf as a given external parameter, we solve the linear equations (5.471 for 
f± and find (5.441 with (5.451; B fol lows as a function of p± by a sin gle iteration. 

We now take 6* < 1. By A± = 0(1 / 6>) wh ile by ( |5.38| w e have Cf = 

0(1/6*2) assuming (/+ + /_)e = 0(1) ^ 0. Thus, (|5.45| leads to (|5.46D. 



Note that the standard Gibbs-Thomson formula (1.3 1 is not applicable here since 
p* — (and hence pq — 0)- However, a linear- in-K term in f± is present, giving rise 
to a "generalized" stiffness 6 that is not bound to a reference density . 

6. Conclusion. The Gibbs-Thomson formula and stiffness of a step edge or 
island boundary were studied systematically from an atomistic, kinetic perspective. 
Our starting point was a kinetic model for out-of-equilibrium processes [71[S]. The 
kinetic effects considered here include diffusion of edge-atoms and convection of kinks 
along step edges, supplemented with mean-field algebraic laws that relate mass fluxes 
to densities. Under the assumption that the model reaches a kinetic steady state 
with straight steps, the step stiffness is determined by perturbing this state for small 
edge curvature and Peclet number P with |k| < 0{P), and applying the quasi-steady 
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approximation . A noteworthy result is that for sufficiently small 6, 0{P^^^) < 6* <C 1, 
the step stiffness behaves as f3 = 0{l/d). This behavior is in qualitative agreement 
with independent calculations based on equilibrium statistical mechanics |37l 139] . 

Our analysis offers the first derivation of the step stiffness, a near-equilibrium 
concept, in the context of nonequilibrium kinetics. The results here are thus a step 
towards a better understanding of how evolution out of equilibrium can be reconciled 
with concepts of equilibrium thermodynamics for crystal surfaces. Furthermore, this 
analysis provides a linkage of microscopic parameters, e.g. atomistic transition rates 
and coordination numbers, to mesoscopic parameters of a BCF-type description. This 
simpler description is often a more attractive alternative for numerical simulations of 
epitaxial growth. 

There are various aspects of the problem that were not addressed in our analysis. 
For instance, it remains an open research direction to compare our predictions with 
results stemming from other kinetic models [H [31 [Hj. The existence of a kinetic 
steady state with straight edges, although expected intuitively for a class of initial 
data, should be tested with numerical computations. Germane is the assumption of 
linear-in-K corrections in expansions for the associated densities. Our perturbation 
analysis is limited by the magnitudes of k and P; specifically, \k\ < 0{P). The 
formal derivations need to be re- worked for k > 0{P) as P ^ 0. The kinetic steady 
state here forms a basis solution for our perturbation theory, and is different from 
an equilibrium state. At equilibrium, detailed balance implies that the fluxes /-|_, /_ 
and each of the physical contributions (terms with different coordination numbers) 



in ( 2.23 H( 2.25 1 for w, g, and h must vanish identically |7]. An analysis based on this 
equilibrium approach and comparisons with the present results are the subjects of 
work in progress. Generally, it also remains a challenge to compare in detail kinetic 
models such as ours with predictions put forth by Kallunki and Krug with regard to 
the Einstein relation for atom migration along a step edge |23j . Our underlying step 
edge model is based on a simple cubic lattice, and it does not include separate rates 
for kink or corner rounding. 

Lastly, we mention two limitations inherent to our model. The mean-field laws for 
the mass fluxes are probably inadequate in physical situations where atom correlations 
are crucial. The study of effects beyond mean field, a compeUing but difficult task, lies 
beyond our present scope. In the same vein, we expect that the effects of elasticity [51 
1251 133] will in principle modify the mesoscopic kinetic rates (attachment-detachment 
and permeability coefficients) and the step stiffness. The inclusion of elastic effects in 
the kinetic model and the study of their implications is a viable direction of near-future 
work. 
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Appendix A. Step edge coordinates and basic relations. In this appendix 
we describe several coordinate systems for an island boundary, thus supplementing 



the formulation of section 2.1 Consider step boundaries that stem from perturbing a 
straight step edge parallel to the a;-axis; see Figure pTl] Three associated coordinates 
and generic densities and longitudinal velocities (along the step edge) are defined as 
follows. 

• Fixed (x-) axis: Variable x, velocity w, density ^ (^ = or fc). 
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• Lagrangian: Variable a, velocity^JI^, density^S. 

• Arc length: Variable s, velocity W , density S. 



The vector- valued normal velocity of the boundary is v n, where n is defined in (2.8 1. 
We note the relations 



(A.l) 



yl = i/cos( 



(A.2) 



= 1 XsXss + VsVss = . 



By use of the Lagrangian coordinate a, we denote 

d 
dt 



(A.3) 



dt\a, 



dt ■■= dt\x , 



i.e., d/dt is the time derivative with the spatial variable a held fixed. Because the 
arc length, s, is only defined up to an arbitrary shift, we choose not to use a time 
derivative with s held fixed; instead, we use d/dt in conjunction with the s derivatives. 
Thus, the interface velocity un in the different coordinates is given by 



(A.4) 



dt\a{x,v) = wn = v{ys, -Xs) , 
dt\x{x, y) ^vh + uiT = (0, U2) 



in which t is defined in (2.8) and 

(A. 5) ui = —vys/xs = —vi&nO , U2 = —v/xs = —v/cosi 

The tangential derivatives and the time derivatives are related by 



(A.6) 



da = Sads , dx = Sxds = (1/ COS 6*) ds , 

dt\x = dt\a + {dt\xa)da = dt\a + {Ul/Sa)da ■ 



We now use these relations to state transformation rules involving the (9, t) vari- 
ables; see section [2T3] for their applications. We assume that 6* is a monotone function 
of the coordinate x and the arc length, s. Useful derivatives in x and t are 



dx ^ Sxds = Sx 9s de 



(A.7) 
(A.8) 

5tU = -7- + {,dt\xa)da 
dt 

where Sx — l/cos6', 9s - 
(A.9) 



cos 9 



d, 



cos 9 cos 9 



de 



-r + idt\xCt)sa ds — y tan 9 9sdg — — 

dt ^ ^ ' dt dt 

—K, and 

d 
dt 



vKta,n9 da 



dt 



dt\e + {dtUe)de , 



(A.IO) 

By ( |A.9[ ) and ( |A.10[ ) we have 
(A.ll) 



dt\a9 = —9 ^ Kve 
dt 



Kve dg 
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Thus, (A. 8) becomes 

(A.12) dt\o. = dt\e + Kivt&ne + ve)de . 
In particular, 

(A. 13) dt\xp{0) ^ k{v tan 6 + ve)dep , 

for any differentiable function p{9) (9t|ep = 0). 

We close this appendix by deriving relations for the densities and velocities along 
a step edge in the different coordinates. If ^, S and S denote line densities of the 
same atom species in x, a and s, we have ^ dx — "Bda — "Bds. Thus, 

S = Xa^ = XsSa C = (cos 9)Sa^ , 

(A.14) E^x4^{cos9)^^E/s^ . 

Next, we derive corresponding relations for the longitudinal velocities w, W and W. 
If the position of a moving point is X(t), S{t) or A(t) in the x, s and a coordinates, 
respectively, then the velocities in these coordinates are related by 

Xt^v/us^v/ sin 9 , 
W ^ At^{w - vys)/xa , 
(A.15) W^Wso,. 

Appendix B. Identities for step edge motion. In this appendix we state and 
prove three propositions pertaining to motion along a step edge. Some of these results 
are used in relation to section |2.3| and in appendix [C] in order to derive alternative 
equations of motion for edge-atom and kink densities. 

Proposition B.l. In the {x,t) variables, the step edge velocity v satisfies 

(B.l) at(tan6')-l-9^(u/cos6') = . 



Proof. We proceed by direct evaluation of the derivatives appearing in (B.l I 
First, we calculate the time derivative in terms of s derivatives via the relation 

(B.2) dtUt'dn9= [dtlo.- —d^]t&n( 

SaXg 



By (A. 2 1 of appendix [A| we evaluate separately each term in the right-hand side: 



9t|atan6' = dt{ya/Xa) = {yatXa - Xatya)/Xa 

[{-VXa/Sa)aXa " {vy^/ Sa)aya]/x'i = -[{vXs)sXs + {vys)sys\/xl 

(B.3) = -[vs{xl + yl) + v{xssXs + yssys)]/xl = -Vs/xl = -vjcos^ 9 , 



2 



da tan 9 = Sads{ys/xs) = Saixsyss - ysXss)/x. 

(B.4) = Sa[Xs{-XsXss/ys) - ysXss\/xl = - '^"'^^ SaXss_ 

VsXg 



Second, we address the spatial derivative in (B.l I: 



B.5 dx{v cos9) = -ds = — ^ H a^s — = — ^ - v — - 

cos 9 cos 9 cos^ cosb* a;, cos^ cos'' ( 
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Equations (B.2|-(B.5l combined yield (B.ll, which completes the proof. 



Proposition B.2. If the line (step-edge) density ^{x,t) satisfies 
(B.6) dt^ + d^{w0^b, 

in the {x,t) coordinates, then the following relations hold in {a,t) and {s,t): 
(B.7) {d/dt)E + da{WE) = B ■.= bxc , 



(B. 



{d/dt)E + ds{WE) + vkE = B := B/s^ 



where dt = dt\x and d/dt = dt\a- 

Proof. We proceed by direct evaluation of the left-hand side of (B.6 1 in the {a,t) 
and {s,t) coordinates. First, we prove (B.7l. By (A. 14) we have ^ = E/xa and 

— {Wxa + vys)E/xa — WE + {vya/xaSa)E. Thus, the time derivative in (B.6 1 
becomes 



i9tUC = dt\a 



(B.9) 



Similarly, the spatial derivative in (B.6) reads 



= a„ [WE + ^E 



(B.IO) 



1 



d<.{WE) 



vya „ 



The combination of iBM and (B.IO I yields 



dt^ + d^iwO - {l/x^ME + d^iWE)] 



(B.ll) 



dt 



=: Ci -y + C2 Va 



X(\Sr 



The last expression is justified by the identity 



dt{l/Xa) = -{Xt)a/xl^ = ~{vys)sSa/xl^ = {-VsysSa " VyssSa)/x'i 



(B.12) 



ya 



yss'^a 



Next, we show that the coefficients Cj in (B.ll) vanish identically. To this end, we 



convert the related s derivatives to a derivatives via the identities 

(B.13) yss = {ya/Sa)a/Sa = yaa/s^ - yaSaa/sl 



'a 5 



(B.14) 



vya 



' + " 2r^ 
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It follows by (B.ll I that Cj are 
(B.15) C2 = 



1 Va 



, 



(B.I6) Ci = -^ 



1 Uaa-^a^a ycti^-^a^aa. H" ^a^^aa) 



= 



which in view of ( B.6 1 and (B.ll I yield (B.7l. 

To derive (B.8I, we first note that WE. — WE.. For arbitrary ai, a-i we consider 
the integral 



Oil 



ax 
a2 



{dt\a^ — ^VKSa) da = / [—da{W^)+B~-^VKSa]da 



[-ds{WE) + B- EvK]sa da= [-~ds{WE) + B- Evk] ds , 



(B.17) 



where we used dt\aSa — vnsa and ( B.7[ ). Equation (B.8I follows directly, thus con- 
cluding the proof. □ 
Proposition B.3. If the line density S,{x,t) satisfies 

(B.18) dt^ - d^idd^O ^ b , 

then the following relations hold: 



(B.19) 
(B.20) 



{d/dt)E - d,{Dd,0 + da{UE) = B := bx^ 



{d/dt)E - ds{DdsE) + ds{UE) + vkE = B:= b/s^ 
where dt — dt\x, d/dt — dt\a, and 
d 



(B.21) 



D := 

x1 



U := dsa 



Xn 



vys 



vys 



Proof. Equations (B.19 1 and (B.20 1 follow directly from Proposition B.2 by setting 
w — —^xl^- Indeed, with this substitution we have 

w - vys 



= ^w- ^vys -dxf - fvys 



— wys = - - 



da 



vys 
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where dgSa — Saa/sa was used in the last hne. Equation (B.19I comes from Propo- 
sition B.2 in view of (B.22). Equation (B.20I stems from Proposition B.2 via (B.23I, 



which completes the proof. □ 

Appendix C. Edge-atom and kink motion in s coordinate. Next, we apply 
the results of appendix [b] to transform evolution laws (2.171 and (2.191 for (j) and k 
to the {s,t) coordinates. 

By use of Proposition B.2 the convection equation (2.19) becomes 

(C.l) {d/dt)K + ds[W{Kr - Ki)] + vkK = 2Fk , 

where 

w — vz/s w — V sin 9 



(C.2) 



Fk (.9 - h)xs ^ (g-h) cose , W 



In the above, K = kcosd, Kr = kr cos 0, Ki = ki cos 6 and W are defined along the 
edge arc length (s) according to the notation of appendix [A| 

Proposition B.3 of appendix [b] converts the diffusion equation (2.17) for (j) to 

(C.3) (d/dt)^ - d.iDEds^) + 9,(C7$) + VK^ = , 

where 

De 



(C.4) 



F. 



f+ + f-- fo cos 



De 

X? 



cos^ 



(C.5) 



U := De Sa 



vy± 



and a is the Lagrangian step coordinate; see appendix [a] Here, <& is the edge-atom 
density defined along the edge arc length. Note that the transformed equation (C.3 1 
contains a drift term, which is absent in (2.17) if x is simply replaced by s. 



Appendix D. Step edge velocity. In this appendix we derive (2.28) in the 
form of a proposition; cf. equation (2.12) in [8 . 

Proposition D.l. The net flux /o oj terrace and edge-atoms to kinks is 



(D.l) 



/o 



:(1 -I- K(j) COS 0) . 



Proof. We apply mass conservation, revisiting the derivation in [T. The starting 
point is the change of the total number of adatoms on a terrace, which is balanced 
by: (i) the step edge motion, (ii) the change of the number of edge-atoms, and (iii) 
the flux rate of deposited atoms. Hence, 



(D.2) 



d 
dt 



pdA 



V ds 



d 
dt 



<^ds-FA 



where A is the area of a single terrace and F is the step boundary. 

Next, we find alternative expressions for the terms ^ / pdA and ^ / ^ds. First 



integration of the diffusion equation (2.13) for p yields 



(D.3) 



dt 



pdA=- {f+ + U)ds + FA 
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Second, direct differentiation of J ^ ds with respect to time gives 

d 



(D.4) 
by using 
(D.5) 



dt 



$(is= / {dt\s^ + Kv^) ds 



A. 
dt 



ds = dt\aSa da = {dt\aSa) da = {KVSa)da = Kv ds 



By combination of (D.2|-(D.4| we obtain 



V ds 



(/+ + /_)ds + 



dt 



* - {dt\aS)ds^ + KV^ 



(D.6) 



(/+ + /_) ds + [F^- ds 
r Jr 



where we invoked the evolution equation (C.3l for $ in {s,t) coordinates and defini- 
tion (C.4| from appendix [C| Thus, via integration by parts, ( |D.6 1 becomes 



J vds ^ - J foXsds 



foXsds+ / $ s^^daidtlas) ds 



(D.7) 

Hence, we have 
(D.8) 



I foXs ds 



<I>K V ds 



which is identified with (2.28) and, thus, concludes the proof. 



Appendix E. First-order perturbation theory. In this appendix we describe 
in the form of a proposition the basic linear perturbation for the equations of motion 
along a step edge. This theory is used in section |5] 

Proposition E.l. Let and k be functions of {0,t) that satisfy 



(E.l) 



Mj{(j), k, K(f>g, Kkg) —Q j = 1, 2 , 



where Mj{(j),k,r],Q are differentiable. If {2.29) holds, where (j)^^) and k^'^^ solve 
(E.2) 



Mj(0("\fc(o\O,O) = O , 



then (p^^) and k^^^ are 

(E.3) 

where 
(E.4) 



X)0 
1^ 



V = 



d^M2 dkM2 
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(E.5) 



{d^Mi)deko 



dkMi 
dkM2 



-{dcMi)d0(j)o-{d^Mi)deko 
-{dQM2)de<f>o-{d^M2)deko 



(E.6) 

and the derivatives oj Mj{<f>,k,'q,C,) are evaluated at A:(°), 0, 0). 

P roof. Equations (E.3|-(E.6) follow directly from the Taylor expansion of for- 
mula (|ki} at (0(°\fc(°^O^), 



(E.7) = d^Mj + {Kk^^^)dkMj + {Kdg(t>^-°^)d^Mj + {Kdgk^^^)d^A'l 



3 ' 



where use was made of dK2l. The 2x2 linear system for fc^^)) leads to |K3|. □ 
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